Plot Utilities

8169356c-68fa-4a08-85ed-6695a9826283 13 15 12 decimal

&[DATE] &[TIME] - &[FILENAME]

Plot Utilities

vx vy char size clr plotG n vx length : plot vx 1 el vy 1 el char size clr augment : i 2 n range plot plot vx i el vy i el char size clr augment stack : for plot 41 line :

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

Bar plot unicolor style, bar width user. data dt BarSize v data 2 col : t v bar t dt - 0 t dt - v t dt + v t dt + 042 mat : j 1 data rows range j 1 ≡ B data j 1 el data j 2 el bar : B B data j 1 el data j 2 el bar stack : if 11 line for B 41 line :

LM(init,Tol,x,y) . . Carlos 2018

init Tol x y LM n β length : m x length : 21 line x β ia x β f : i 1 m range ΦΦ i el x i el β ia y i el - : for β Φ ΦΦ : 21 line i 1 m range j 1 n range β j el β j el h + : jjh1 x i el β ia y i el - : β j el β j el 2 h * - : jjh2 x i el β ia y i el - : jjh0 i j el jjh1 jjh2 - 2 h * / : β j el β j el h + : 61 line for 11 line for β derfhφ jjh0 : 51 line cont 0 : dnor 1 : update 1 : 31 line h Tol : μ 0.01 : μ1 100 : 31 line 12 mat dnor Tol > update 1 ≡ f init Φ : J init derfhφ : jtj J transpose J * : res f norme (2^: 41 line if 11 line a jtj μ n identity * + : s1 a 1 -^J transpose f * (* : b1 init s1 - (: f1 b1 Φ : SSD f1 norme (2^μ1 / : 51 line SSD res < init b1 : f f1 : 12 mat dnor s1 normi init normi / : cont cont 1 + : Sol cont 1 el cont : z 1 n range Sol cont z 1 + el init z el : 11 line for Sol cont n 2 + el SSD μ1 * : Sol cont n 3 + el μ : Sol cont n 4 + el dnor : μ μ 10 / : update 1 : 101 line μ μ 10 * : update 0 : 21 line if 31 line while β init : 41 line :

# exp

2. Alternate global fit ... SS 6179

x 0 0.51 0.82 1.02 1.23 1.43 1.85 2.05 2.36 4.52 9.22 111 mat :

y 0.78 0.8 0.92 1.14 1.26 1.28 1.22 1.16 1.12 1.08 1.08 111 mat :

Parmacokinetics model function ... Mathcad fails !!! x β f β 2 el β 1 el β 2 el - 1 x β 3 el / (β 4 el^+ / + β 1 el β 2 el - 1 x β 3 el / (β 4 el^+ (2^/ x β 3 el / (β 4 el^* β 4 el * - :

β β 1 el β 2 el β 3 el β 4 el 14 mat transpose :

... Prepare L-M algo ... components fetched by algo

Very robust wrt initials init 0.5 11114 mat transpose :

Typical 'Tol' this model Tol 0.0001 :

t0 1 time :

=================================

β init Tol x y LM :

recover @ 'Tol' β β 4 round vectorize :

1 time t0 - 23 s *

β 0.7773 1.0727 1.2781 4.8008 41 mat

0.7773 1.0727 1.2781 4.8009 41 mat

resd i 1 x rows range Δy i el y i el x i el β f - eval : for x Δy 2 round vectorize augment 21 line :

UV x y augment : data UV o 5 black plot : x c 10 x ≤ (x 10 ≤ (& cases : Dev resd 0.1 BarSize : 41 sys

Superb L-M fit to Pharmacokinetics ... undisclosed from Art.

1.08 x β f x c * data 3 1 sys

resd 00 0.51 0 0.82 0.02 - 1.02 0.02 1.23 0 1.43 0.01 - 1.85 0.01 2.05 0.01 - 2.36 0.01 - 4.52 0 9.22 0.01 112 mat

0 Dev 2 1 sys

Decimate overpopulated data set

data n Decimate collect every 'n' pairs range 1 data rows n range : k 1 : i range u k el data i 1 el : v k el data i 2 el : k k 1 + : 31 line for u v augment 51 line :

D2 ... D3

D2 00 0.021 2.22 0.04 3.375 0.059 3.789 0.079 3.782 0.099 3.538 0.12 3.164 0.141 2.754 0.164 2.322 0.188 1.919 0.213 1.552 0.237 1.253 0.263 0.991 0.29 0.77 0.318 0.595 0.346 0.456 0.377 0.341 0.406 0.257 0.437 0.192 0.47 0.14 0.509 0.096 0.545 0.068 0.591 0.044 0.645 0.026 0.918 0.002 1.125 0 1.154 0272 mat :

D3 0 10.861 0.021 6.056 0.04 3.327 0.059 1.687 0.079 0.706 0.099 0.133 0.12 0.187 - 0.141 0.338 - 0.164 0.395 - 0.188 0.394 - 0.213 0.362 - 0.237 0.32 - 0.263 0.272 - 0.29 0.226 - 0.318 0.185 - 0.346 0.149 - 0.377 0.118 - 0.406 0.093 - 0.437 0.073 - 0.47 0.056 - 0.509 0.041 - 0.545 0.031 - 0.591 0.022 - 0.645 0.014 - 0.918 0.002 - 1.125 0 1.154 0272 mat :

Error stem plot for QuickPlot unicolor data Error stemX 1 data rows range : i 1 stemX rows range stem i el data i 1 el data i 2 el data i 1 el 022 mat : for i 1 stem rows 1 - (range stem i el stem i el stem i 1 + (1 el 1 el 012 mat stack : for V stem 1 el : j 2 stem rows 1 - (range V V stem j el stack : for V V data data rows row stack : 71 line :

XY D2 :

raw XY 3 Decimate o 5 black plot :

raw